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■ We calculate relativistic Fermi liquid parameters (RFLPs) for the description of the 

(N 

p) properties of dense nuclear matter (DNM) using Effective Chiral Model. Analytical 
D ■ 

^/^ I expressions of Fermi liquid parameters (FLPs) are presented both for the direct and 

■ exchange contributions. We present a comparative study of perturbative calculation 
with mean field (MF) results. Moreover we go beyond the MF so as to estimate 

Ph' the pionic contribution to the FLPs. Finally, we use these parameters to estimate 

Qj • some of the bulk quantities like incompressibility, sound velocity, symmetry energy 

^ : 

etc. for DNM interacting via exchange of a, uj and vr meson. In addition, we also 
, calculate the energy densities and the binding energy curve for the nuclear matter. 

I Results for the latter have been found to be consistent with two loop calculations 

^Sj ■ reported recently within the same model. 
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I. INTRODUCTION 

One of the most exciting field of contemporary nuclear research has been the studies 
of the properties of dense nuclear matter (DNM). Such studies are important both in the 
context of laboratory experiments and nuclear astrophysics. Therefore, several attempts 
have been made in recent years to ascertain the properties of nuclear system at densities 
higher than the normal nuclear matter densities [l, 2, sl, ^, I^, B]- 

The suitable description of nuclear matter at such high densities is provided by Quantum 



* Electronic address: | kausik. pal@saha.ac.in| 



Hadrodynamics (QHD) 



ll| to study 



Historically, QHD was developed by Walecka 
the properties of neutron star where the nucleons are assumed to interact via the exchange 
of a and uo mesons. In this model starting with interacting Lagrangian the relativistic field 



equations are solved by making MF approximation where the meson fields are replaced 
by their vacuum expectation values. Subsequently, starting from the same model Chin 
developed a full diagrammatic scheme and showed that the MF results can be obtained 
by making Hartree approximation i.e. by retaining only the direct terms in a relativistic 
field theoretic approach [l^ • In the same work, it was also shown how exchange corrections 
can be made and analytical expressions can be found for the energy density and related 
quantities by making some long range approximation for the n — n interaction. Since then 



the QHD has undergone a series o 



reader to ref. 
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developments which we do not discuss here but refer the 



19|. 



The most recent model which we use here for the des crip tion of dense nuclear system is 
provided by the Chiral Effective Field theory (chEFT) 20|, l2l|]. It might be recalled here, 
that, in such a framework, the explicit calculation of the Dirac vacuum is not required, 
rather, on the contrary, here, the short distance dynamics are absorbed into the parameters 
of the theory adjusted phenomenologically by fitting empirical data. For detail discussion 



refer the reader to ref. 
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24( 1 . Recently this model has been applied 5 



to calculate 



the exchange corrections by evaluating nucleon loops involving a, uj and tt as intermediate 
states, which we address here partly. 



Our approach here is to study the dense nuclear system in terms of relativistic 



ermi 



liquid parameters (RFLPs). Such an extention of the Fermi Liquid theory 25|, |26|| was 
first made by Baym and Chin in ref.j23|. It should, however, be noted that the calculation 
presented in ref. 27|] were performed perturbatively where the original QHD model was used. 
It should, however, be noted that the first application of Fermi liquid theory (FLT) to study 



the nuclear system was due to Migdal 
unbound nuclear matter and finite nuclei 
the nuclear shell model 



28]who used FLT to investigate the properties of 
29!]. FLT also provides theoretical foundation for 



291] as well as nuclear dynamics of low energy excitations 26|, |30| . The 



connection between Landau, Brueckner-Bethe and Migdal theories was discussed in ref. 31]. 
While these are all non-relativistic calculations, the relativistic calculations involving RFLT 
are rather limited. 

After the original work of the relativistic problem was revisited in 32] where one 
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starts from the expression of energy density in presence of scalar and vector meson MF and 
takes functional derivatives to determine the FLPs. The results are found to be qualitatively 



different than the perturbative results 27|, |32|. Moreover, besides a and uj meson, ref. 32| 



also includes the p and vr meson and the model adopted was originally proposed by Serot 
that incorporates pion into the Walecka model. The latter, however, do not contribute to the 



parameters presented in 



32| as the calculation was restricted only to the MF level where pion 



fails to contribute. On the other hand, the Mi 



models of the nuclear force calculated in ref. 



gda^ 



33 



parameters using one-boson-exchange 



34| . in which a comparison of relativistic 



and non-relativistic results have also been studied. 

In the present work, we use a model, where we have pions and we extend the calculation 
beyond MF to include the pionic contributions into the FLPs. Furthermore, we evaluate 
and compare the perturbative results with MF approximated results within the framework 
of the present model. In addition we also calculate various physical quantities like incom- 
pressibility, sound velocity and symmetry energy etc. Moreover, the results are compared 
whenever possible with the previous calculations by taking suitable limits. For instance, 
the exchange energy, we compare results calculated within the present scheme with a more 
direct evaluation of the loop diagrams like in ref. 231 ]. 

This paper is organized as follows. In Sec. II, we will depict brief outline of the formal- 
ism of FLT. We find the analytic expressions for the FLPs both for direct and exchange 
contributions in Sec. III. Subsequently, we determine chemical potential, energy density and 
various other thermodynamic quantities like incompressibility and sound velocity. Sec. IV, 
is devoted to calculate isovector LPs to which involves the vr meson contribution, and used 
to express the symmetry energy. 



II. FORMALISM 



In FLT total energy E of an interacting system is the functional of occupation number 
Up of the quasi-particle states of momentum p. The excitation of the system is equivalent 
to the change of o ccup ation number by an amount 5np. The corresponding energy of the 
system is given by [26|, 1271] , 
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(27r)3 (^27r)3-^^'''^'''^ 



(1) 



where is the ground state energy and s is the spin index, and the quasi-particle energy 
can be written as, 



e =6° +V 



d^p' 



(27r 



(2) 



where 5°^ is the non-interacting single particle energy. The interaction between quasi- 



particles is given by fps,p's', which is defined to be the second derivative of the energy 
functional with respect to occupation functions. 



6^E 



ps,p's' — e- t- 



Since quasi-particles are well defined only near the Fermi surface, one assumes 



(3) 



fl + Vf{p-pf) 



and p ~ p' ^ pf. 



Then LPs //s are defined by the Legendre expansion of fps,p's' as |26|, |27|, 



(4) 



fi 



21 + 1 



(5) 



where 9 is the angle between p and p', both taken to be on the Fermi surface, and the 
integration is over all directions of p 27|]. We restrict ourselves for / < 1 i.e. /o and /i, since 
higher / contribution decreases rapidly. 



Now the Landau Fermi 
scattering amplitude via 
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iquid interaction fps,p's' is related to the two particle forward 
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FIG. 1: Schematic diagrams of direct contribution to forward scattering amplitude. Nucleons are 
represented by solid lines, cr, uj and tt mesons are denoted by dotted, wavy and dashed lines 
respectively. 

where M is the mass of the nucleon and the Lorentz invariant matrix M.ps^pis' consists of 
the usual direct and exchange amphtude, which may be evaluated directly from the relevant 
diagrams as shown in Fig. [1] and Fig. [2l 
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FIG. 2: Schematic diagrams of exchange contribution to forward scattering amplitude. Nucleons 
are represented by solid lines, cr, uj and vr mesons are denoted by dotted, wavy and dashed lines 
respectively. 

The spin averaged scattering amplitude (/pp') is given by {2?! . 



4 ^ e^, e;, 



(7) 



The dimensionless LPs are Fi = N{0)fi, where A^(0) is the density of states at the Fermi 
surface defined as, 
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9s9ipj f dp 



~ 27r2 ■ 

Here gs , gi are the spin and isospin degeneracy factor respectively. 

In the above expression {dp/dep)p=p^ is the inverse Fermi velocity {vj'^) related to the 
FL parameter Fi, 

vj' = idp/dep)p=p^. = iii/pf)il + Fi/3). (9) 



With the help of Eq.® and Eq.(Q one writes |32| 



ej=fx{l + ^F,). (10) 

To compare Eq.(l9]) and Eq.( |TOl) with the well known non-relativistic expressions one has 
to put Ef = M* and fi = M. 



III. CHIRAL LAGRANGIAN AND LANDAU PARAMETERS 



We adopt the non-linear chiral model to calculate the FLPs and consequently estimate 
various quantities of physical interest like effective chemical potential, sound velocity, in- 
compressibility, symmetry energy etc. Here all the fields are treated relativistically 27|. By 



retaining only the lowest order terms in the pion fields, one obtains the following Lagrangian 



from the chirally invariant Lagrangian 
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23| 



1 



(11) 
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where cu^jy = S^cu^ — dyU^, n= i(r ■ $7r) and r is the isospin index. Here \1/ is the nucleon 
field and and are the meson fields (isoscalar-scalar, isoscalar- vector and isovector- 

pseudoscalar respectively). The terms 6C and Cnl contain the non-linear and counterterms 
respectively (for explicit expression seej23]). Note that in this work, the convention of [l^ 
is used. 



Pert ur bat ive calculation 



Let us calculate the LPs perturbatively due to the exchange of scalar and vector mesons 
between the nucleons The direct contribution [see Fig. [1]) is given by 27] 



rdir,a 
J pp' 

rdir,uj 
J pp' 



I P.P' 



(12) 



where = a/p^ + M^. Now with the help of Eq.([5]) and Eq. (fT2|) . the LPs become 



f 
Jo 



dir,(7 



Jo 



dir,uj 



9cu 



(13) 



and 



One may neglect t 
developed by Matsui 



rdir,(j 

Ji — u 

pdir,uj _ P} 
•'1 rnll^' 



le contribution of f'^^^''^ as discussed in ref. [2[. A better approach was 



(14) 



.a. 



321 ] where the magnetic interaction is included which reduces the value 



of/i 



dir,L 



One may now, for the direct contribution plug in f^l?'^ and /^J''^ in Eq. ([T]) and Eq. ([2]) to 



pp 



edir,L 
I pp' 



obtain the energy density and the SPE spectrum, respectively. The SPE spectrum is given 

by y 



^d^r ^ ,0 9l 



(15) 
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Here p and Us are the baryon and scalar density given by 



P = 9sgi 



and 



(16) 



9sgi 



M 
4v]-2 



The energy density for direct contribution is 27 1 



(17) 



dir rpQ . ^ iJu 2 ^ Ua 2 



The chemical potential is 



2 m 



2 ml 



(18) 



dir- 



dE 



'dir 



dp 

JP J^s^r- 

op 

^gl 9lM 

ef + —p o—^s- 

m^, mi Ef 

Lu U J 



One can derive the same result directly from Eq. (ITSl) as p = Ep 



p=Pf 



(19) 



B. FLPs in mean field model 



It is well known that in the MF approximation, one replaces the mesonic fields by their 
vacuum expectation values viz. a ^< a >= ao, uj ^< uj >= S^qU^. The pion, however, 



fails to contribute at the M 



density can be written as lOj 



level as < TT >= 0. In the MF approximation the energy 



E 



MFT 



2m,^ 



2 m 



(20) 



In the above equation M* denotes the effective nucleon mass to be determined self con- 



sistently 
form 



23|, 



32] • With the help of Eq.([n]), the interaction parameter takes the following 
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rMFT 
J pp' 



„2 Tur*2 



Here = a/p^ + M*^ and 



(21) 



C(M* 




d .j*2 9sgiv f 

^2 27r2 



4vt 



In 



1 + ^/ 



(22) 

where f / = Pf/ipj + M*^)^/^, is the relativistic Fermi velocity. The inverse part of Eq.( l2T| 
reduces the magnitude of interaction parameter compared to what is obtained in absence of 
the MF Eq.([T2]) 



32 



34| 



The LPs as defined in ref. 



32| are 



f 
Jo 



MFT,a- 



^ T~ 

mi ej 



f 

Jo 



MFT,LO 



1 + %E,-^ 



*(p|+M*2)3/2 



(23) 



When we evaluated the above Eg. (1231) . we neglect the "magnetic interaction" between the 
quasiparticles which is induced by the microscopic currents. In presence of current density 
32| 



/i 



MFT,uj 



ml ej 



1 + ^E-^ 



M 



*2 



{pf + M*2)3/2 



(24) 



Clearly, the current contribution reduces the value of f^^^'^'^ ^ Previously we showed in 
Eq. flT^ and Eq. f[TH]l the SPE spectrum and energy density in absence of MF, but in presence 
of MF, SPE is given by 
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32| 



.MFT 



9uj I / 



M 



*2 



(25) 



Therefore, 
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In the low density limit, Eq.fl2Sl) reduces to Eq.flT^ as M* = M — -^rj^i^ 



M 



i"'*(p2+Af2)l/2 



. It 



is to be noted that in the MF approximation scalar meson contribution is absorbed in the 
effective mass does not appear explicitly as in Eq. (fT5|) . Another interesting difference is also 
noticed in the expressions for the total e nerg y densities given by Egs. ffTSj) a nd fl20|) . Note 
that, our MF result is consistent with ref. 32(] but differs with that of ref.[ll|. 
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FIG. 3: Chemical potential for direct contribution with a and to meson exchange in symmetric 
nuclear matter. The dashed and solid curve represent the perturbative and MF results, respectively. 
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FIG. 4: Total energy from direct contribution with a and co meson exchange in symmetric nuclear 
matter. The dashed and solid curve represent the perturbative and MF results, respectively. 



In Fig(l3]) we present the comparative study of the chemical potential obtained pertur- 
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batively and with MF calculation. At low density they tend to merge, while at higher 
density MF results differ significantly from the perturbative result. Numerically /i*^"^^ 
and /i^'^'' are given by 861.07 MeV and 832.64 MeV respectively at normal matter den- 
sity (po = 0.148fm~'^). For the numerical estimate we adopt the coupling parameter set as 
designated by MO A in ref. 23|]. In Fig(jl]) we compare the results for total energy obtained 
from perturbative and MF calculation. This also shows at low density they tend to merge, 
while at higher density MF results become larger than the perturbative results. This is 
easily understood from Eq. flTH]) and fl2Ul) . At saturation density numerical values are given 
by 886.43 MeV and 893.31 MeV for perturbative and MF calculation respectively. 

Now we consider the exchange modification over the MF. Evaluating the exchange dia- 
grams (Fig. [2]), we obtain the interaction parameter as {27] 



rex, a 
J pp' 

rex. 
J riY}' 



P.P'+M* 



ex,LU 
pp' 



■4ipn;(p_p,)2_„ 
gg, P.P'~2M* 

■^i0^(p_p/)2„„ 



With the help of Eq.([5]), the LPs for scalar meson exchange reads as 



and 



rex, a iter 

Jo — 



o2 /•! p){l-cose)+2M*'^ ^ 
' — d(cosi 



9l 



cos I 



1 2p2(l 

ml - 4M*2 



4p2 



/ 



In 1 + 



4p} 



- rec 



cost 



2M 



*2 



8ej 2p}{l - cos^) +m2 



8e} 



(cos^)d(cos I 



ml - AM 




m„ 



2p} 



4p) 



In 1 + 




Using Eqs.([28]) and ([29]) we have 



(27) 



(28) 



(29) 



rex,a 

Jo 



■ rex,a 



cos^)d(cos^) 




(30) 
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It is this combination i.e. /q — which appears in the calculation of chemical potential 
and other relevant quantities. For massless scalar meson interaction, the above Eq. (l30l) turns 
out to be finite, 



Note that in the limit mo- — > both /o and /i are individually diverge because of the 
presence of (1 — cos^) term in the denominator. In the massless limit such divergences are 
also contained in Eq.f l28l) and Eq.( l29i) . 

The dimensionless LPs Fq and Fi are defined as Fq = A^(0)/o and Fi = N{0)fi, where 
A^(0) is the density of states at the Fermi surface defined in Eq.Q. 
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4p2 



(32) 



and 



9s9i 



ml - 4M*2 





(33) 



Similarly, for vector meson exchange we have 



cos( 



M 



*2 



1 2p}(l 



cos I 



-d(cosi 



-2 + i!^i±|^ln 



(34) 



and 



r/i 




Using Eq. 



16£:jPj 
and Eq. fl35l) we have 



(35) 
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TABLE I: Parameter sets used in this work. 



Meson 


Mass 


Coupling 


a 


m^/M=0.54 


c,^/47r=0.7936 




m^/M=0.8328 


5^/47r=0.9681 



Jo 



As} 



As} 



cos( 



M 



*2 



cos t/j 



mf: 



4p2 



[1 — cos6')d(cos( 



In the limit niij —>■ the above Eq. fl36l) turns into 



(36) 



^ 3 



4p5 



2M 



*2 



(37) 



This expression agrees with the previous calculation by Baym and Chin [27| who arrived 
at this result by direct evaluation of the integral by putting m^^ = in Eq. (!36|) . Here also to 
be noted, in the limit both /o and /i are individually divergent, but the combination 

/o — |/i is finite as observed in the case of scalar (a) meson exchange. 

The dimensionless LPs for vector meson exchange reads: 



and 



T-iex,uj 
^0 



-9sgi 



2M 



*2> 



2p) 



In 1 



4p2 



m 



(38) 



3^' 



-9s9r 



2M 



*2\ 



32n^Pfef 



2pj I \ 



(39) 



In Fig. ([5]) and Fig. ([6]) we present F^^, F^^, Fq°* and F^^ as a function of baryon density 
for symmetric nuclear matter due to a and u meson interaction respectively. It is to be 
noted, that F^ ^^ and F^^ contribute in opposite sign for both a and u meson exchange as 
it is seen from Table ([11]) • We quote few numerical values of F'^ and F'^ in Table ([11]) at 
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FIG. 5: Dimensionless LPs in symmetric nuclear matter for a meson exchange in relativistic theory. 
Fq ^^ , Fq^, Fq°* and Ff^ are denoted by dot-dashed, dashed, soUd and dotted hne respectively. 




FIG. 6: Dimensionless LPs for symmetric nuclear matter for uj meson exchange in relativistic theory. 
Fq^^ , Fq^, Fq°* and Ff^ are denoted by dashed, dot-dashed, solid and dotted line respectively. 



normal matter density (po = 0.148fm ^). It is to be noted that the numerical estimation of 



Fa with MF in our case, differs from ref. 



32l |. This is due to different coupling parameters 



in these two models. 

We now proceed to calculate chemical potential due to the exchange terms denoted by 
/i*^^. As in ref. 27| we have 
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TABLE II: Dimensionless LPs for a and uj exchange at p = pQ. 



Meson 








pex 


a 

UJ 


-8.65 
7.35 


3.61 
-1.91 


-5.04 
5.44 


0.875 
-0.932 



dp 



2%^ fde. 



9s9iPf V dp , 

77/1 + hi 

l^9sgiPf 3 



/o 



and 



(40) 



dep 
dp 



p=Pf 



Pi 9s9iPf fi 



27r2 3 



Now from Eg. (1401) one gets 



(41) 



/id/i 



/iPf 1 
Pf + 9s9i^{fo - g/l) 



dpf. 



(42) 



To calculate fi, it is sufficient to let = £f in the right hand side of Eq. (H2|) . With the 
constant of integration adjusted so that at high density pj ~ £/, Eq. (l42|) upon integration 
together with Eq. (l30l) yield 
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-2y.(4-y2)3/2tan-i 



- yl 



^ yl(A-yl)VTT^ ^^ / ^4^ 



+ 2{yl - Qyl + 6) ln(x + VTT^) 
where x = pf/M* and y^j = m^/M*. 



(43) 



For massless meson limit i.e. at — implies ^ we have 
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FIG. 7: Density dependent of exchange chemical potential in symmetric nuclear matter. a,uj and 
TT mesons are denoted by solid, dotted and dashed line respectively. 



Similarly for vector meson interaction, using Eq. (l36l) and Eg. (1421) one obtains 
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yujVT 

4x 



In ( 1 + — I - AxVl + x2 



+ 2(6 - yt) ln(x + Vl + x^] 



(45) 



Here y^^ = m^/M*. For massless limit of vector meson the expression for chemical 
potential reads as 



— + 9s9i- 



f 



X 



Vl+x^ - 3 ln(x + Vl 



x^ 



(46) 



In low density limit (M* M) for the massless meson exchange we reproduce the 



expression derived earlier 



Once the fi'^^ is determined, one can readily calcu! 
density. For scalar meson interaction it is given by 12|, 



ate its contribution to the energy 
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FIG. 8: Comparison of mean field and exchange results of chemical potential in symmetric nuclear 
matter. Direct contributions are plotted by solid curve and exchange contributions by dashed 
curve. 
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FIG. 9: Comparison of chemical potential with MFT and HF case in symmetric nuclear matter. 
MF and HF results denoted by dashed and solid line respectively. 
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(47) 



where 
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Similarly, for vector meson exchange we obtain 



2567r4 
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In 1 



vl 



— X 



[{x\J\ + x2 - ln(x + a/1 + x^)}^ - x^j + 4 



(49) 



where 



2y^(l/^-2£-8) 



tan 



x^ 



dx. 



x^ 



(50) 



Thus for the case of massless mesons 
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27(1, we have 



-[x^ - -{xVl+l^ - ln(x + vTT^)y 



(51) 



and 



i?-U^_o = --^[x' - -{xv/TT^- ln(x + Vl + x^)Y]. (52) 

Due to presence of pion fields in the chiral Lagrangian we have component in the inter- 
action which acts on the isospin fluctuation. One can derive the quasiparticle interaction 
with isospin dependency by the same procedure as for a and uj meson. Pion, being a pseu- 
doscalar, fails to contribute at the MF level forcing us to go beyond MFT so as to include 
pionic contribution to the FLPs. It is to be noted that, in exchange diagram pion have both 
isoscalar and isovector contribution to FLPs. Detailed calculation for isoscalar contribution 
to FLPs is similar as a and uj meson. For brevity, we present only dimensionless FLPs and 
their contribution to energy density. We also quote their numerical values. 

The dimensionless LPs due to vr exchange are 



19 



0.8 




FIG. 10: Dimensionless LPs in symmetric nuclear matter for pion exchange in relativistic theory. 
Solid and dashed line represent Fq and respectively. 
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USTT^PPfSf 

In the non-relativistic limit ej —>■ M*, one obtains the same expression of F[ as reported 



m 



. In Fig. lfTOll we show the density dependence of Fq and Fi due to pionic interaction. 
Numerically at nuclear saturation density (po = 0.148fm~'^), Fq^'^ = 0.68 and F^^'^ = —0.2. 

In Fig. (fTT]) and Fig. (fT2l) we plot separate and total contribution of Fq and Fi due to cr, u 
and TT exchange respectively. Interestingly, individual contribution to LPs of a and u meson 
are large while sum of their contribution to -Fq°* is small due to the sensitive cancellation of 
Fq and Fq as can be seen from Fig. (fTT]) . Such a cancellation is responsible for the nuclear 



saturation dynamics 
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34| . Numerically, Fq~^^ is approximately 3/2 times smaller than 



Fq as can be seen from Table (lllip . 

The exchange energy density is given by 
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FIG. 11: Dimensionless relativistic LP Fq in symmetric nuclear matter, a, to, vr and total contri- 
bution are denoted by dashed, dot-dashed, dotted and solid line respectively. 




FIG. 12: Dimensionless relativistic LP Fi in symmetric nuclear matter, a, to, vr and total contri- 
bution are denoted by dashed, dot-dashed, dotted and solid line respectively. 
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where = rrij^/M* and 
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FIG. 13: Comparison of separate and total exchange energy obtained from FLT in symmetric 
nuclear matter. cr, w, tt and total contribution are denoted by dashed, dot-dashed, dotted and solid 
line respectively. 




FIG. 14: Comparison of mean field energy and exchange energy obtained from FLT in symmetric 
nuclear matter. E^^^, E^^^, E^^, E'^ and E^^ are denoted by solid, dashed, dot-dashed, dash- 
dashed and dotted line respectively. 



/. = -2^3^!^ / ^ tan- ( ) do;. (56) 



For the massless pion this reads as 
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TABLE III: Dimensionless Landau parameters and chemical potential at p = pQ. Note that, Fq, 
Fi are the dimensionless isoscalar LPs. 



Meson 


Fo 






a 


-5.04 


0.875 


73L89 


LO 


5.44 


-0.93 


501.82 


TT 


0.68 


-0.20 


609.88 



TABLE IV: Chemical potential in MeV from FLT at p = po- 
Meson p^^^ 
a 73L89 
cu 50L82 
TT 609.88 
a + LO + TT 675.63 86L07 952.73 



El 



^^^^*'[xr/-ln(x + r/)]2. (57) 



In Fig. f|T3l) and Fig. ffT^ we show the density dependence of energy due a, uj and tt meson 
exchanges. Numerical values are quoted in Table (jVT). 

It might be mentioned here that in the massless meson limit, Eg. (1511) . (l52l) and ( 1571) can 



be evaluated analytically from two loop ring diagrams of ref. [23|] using Eqs.(54), (55) and 
(56). We have checked and the expression for the energies are found to be consistent with 
each other. With massive meson results are compared numerically. Numerical estimation of 
exchange energy from loop diagram and RFLT are found to be few percent limit. 

Finally we reproduce saturation property of nuclear matter i.e. E/p — M = —16.12 MeV 
at pf = 1.3fm~^ with those energy calculated from RFLPs. 
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TABLE V: MF and Exchange energy in MeV from FLT at p = pQ. 



Meson 






a 


193.86 


40.48 


u 


138.39 


-23.41 






12.49 




FIG. 15: Binding energy graph from FLT for symmetric nuclear matter. 



C. Incompressibility and First Sound Velocity 



In nuclear matter several important relationships exist between nuclear observables and 
the FLPs. The thermodynamical parameters can be expressed in terms of few LPs. For 
example we present the incompressibility [K) and first sound velocity (ci)}?. Is^. 



Incompressibility of the Fermi liquid may be derived as in the non-relativistic theory by 
the second derivative of energy density {E) with respect to the number density p 0, Is^ ]: 



If energy density E is given in terms of number density p, then the expression for incom- 
pressibility or compression modulus in terms of LPs is given by, 
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FIG. 16: Compressibility K in symmetric nuclear matter. 



K 



+ Fo). 



(59) 



Now consider the effect of quasiparticle collision on the collective modes of a neutral 
Fermi liquid. Suppose the frequency of the mode is u, while the quasiparticle collision 
frequency is u. For the limit u << u, many quasiparticle collision takes place during time 
interval uj~^. Then the region is collision-dominated, or hydrodynamic regime 32|]. Under 
this circumstances, organized density fluctuation is possible and hydrodynamic or flrst sound 
waves are generated. Hydrodynamic sound propagates only when the system attains the 
local thermodynamic equilibrium in a time much shorter than the time interval of the sound 
oscillation. 



The flrst sound velocity is given by 27| 



dP 
'dE 



With the help of Eq.(^, Eq.^U^) yields 



dP dfi dp 
dfi dp dE 
p dp 
pdp' 



(60) 



IPf I + Fq 
3/i2l + |Fi 

-.2 



(61) 
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FIG. 17: First sound velocity ci in symmetric nuclear matter. 



Corresponding values of the incompressibility and the first sound velocity are plotted in 
Fig. f|T6l) and Fig. f|T7|) separately with a + uj and a + uj + 71 contribution. It is observed that 
for combined a and uj meson at p < 0.75po; Fq < —1 and the resulting compressibility turns 
out to be negative. While for vr meson together with a and ut meson the same conclusion 
can be drawn at p < 0.55po • This is the region where the attractive interaction due to 
the exchange of scalar mesons overwhelms the repulsive force coming from vector meson 



exchange, and consequently the system becomes unstable [32|. On the other hand, as the 
density increases, the attractive scalar meson exchange force tends to be suppressed by 
the relativistic effect and the net quasiparticle interaction become repulsive. At nuclear 
saturation density (po = 0.148fm~^) we have K = 476.04 MeV and 705.84 MeV for combined 
a + UJ and a + u + ir mesons respectively. The small effective mass is responsible for large 
incompressibility. The first sound velocity Ci = 0.19 for a + u and 0.23 for a + + vr at 
normal nuclear matter density and at all densities Ci is smaller than the velocity of light, 
which is consistent with causality 32 1. 



IV. ISOVECTOR LANDAU PARAMETER AND SYMMETRY ENERGY 

In this section we proceed to calculate isovector LPs due to pion exchange. The isovector 
contribution to interaction parameter is 
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f 

J pp' 



lg\M 



*2 



COS I 



p^p'=p, 2 4/2£2 I 2p2 (1 - COS 9) + m2 
93Mey and = 139 MeV 



(62) 



where g\ = 1.5876 , /„ 
in isovector channel j^. 

Using Eq.Q and (l62l) we derive isovector LPs /g and f[, 



231]. Here —1/2 is isospin factor 



and 



/o 



(1 



cos( 



1 2p^(l — cos 61) + ml 



d(cos I 



*2 



9\M 



-2 + 



In 1 + 



COS^(l 



COS I 



16/2e2 (1 - cos 6) + m2 



d(cos I 



64/2^^2 



2p] 



1 In 1 + 



4p' 



m: 



(63) 



(64) 



The dimensionless LPs Fq = A^(0)/g and Fl = N{0)f[ are given below. Here, A^(0) is the 
density of states at the Fermi surface defined in Eq.®. 



F' 



QsQi 



647r2/.%/ 



ml , / 4p2. 
2pj \ m^ 



(65) 



and 



3 ^ 



9s9i 



9WnM*^ 



2p} 



1 In 1 



4p' 



(66) 



Knowing the isovector LPs, to which here only the pion contributes, one can calculate nu- 
clear symmetry energy. The symmetry energy is defined as the difference of energy between 
the neutron matter and symmetric nuclear matter which we denote as /3 32 1. 

Analytically, the symmetry energy is defined from the expansion of the energy per nucleon 
E'{p, a) in terms of the asymmetry parameter a defined as 36] 
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FIG. 18: Density dependent of isovector LPs in symmetric nuclear matter. 



a = 



P3 

P 



Pn - Pp 
P 



N-Z 
A ■ 



(67) 



We have 



E'ip, a) ^ = E\p) + K^Mc^' + 0{a') + (68) 

P 



and so, in general, 



E' ^3 

sym — y 



1 d^E'jp, a) 

2 dc^ 
1 d^E 



a=0 



2^ dpi 



ps=0 



(69) 



In terms of LPs, the symmetry energy can be expressed as 



/3 



(70) 



where Fq is the isospin dependent LP. Numerically at saturation density (p = po) we 
obtain (3 = 14.57 MeV. So relatively small contribution to /5 comes from one-pion exchange 
diagram [37|. 
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Similar to the hydro dynamic sound corresponding to the total baryon density oscillation, 
we consider hydrodynamic isospin sound which accompany the out-of-phase oscillations of 
proton and neutron density. Isospin sound velocity c'^ is given by 32| 



(71) 



Numerically at saturation density (p = po) "we obtain c'l = 0.17. 



V. DISCUSSION 

In this work we determine the RFLPs using effective chiral model to study the properties 
of DNM. We compare the MF results with the corresponding perturbative results for the 
FLPs. It is seen at low densities, they converge but perturbative results differs significantly 
with the mean field results at higher density. 

We have estimated the pionic contribution to the FLPs which can contribute only in 
exchange diagrams. It is seen that numerically pionic contribution is smaller than the 
corresponding a and u meson results, however, we show this is still important in various 
physical quantities as a and u meson contribute in opposite sign and pionic contribution 
is comparable to the sum of the contributions coming from the scalar and vector meson 
exchanges. Among the other physical quantities we also evaluate the sound velocity and 
incompressibility and symmetry energy etc. The results for the sound velocities gives the 
correct causal limit at extremely high densities. It is also shown that, the system is stable 
with respect to density fluctuation at and above the nuclear saturation density and exhibit 
instability in relatively low density regimes. 

Here we also have calculated exchange energy densities by using FLPs calculated relativis- 
tically. The results are found to be consistent with what are obtains by directly evaluating 



the nucleon loops [23|, this gives further support to the present approach. 
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